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Abstract 

We have studied the Cabibbo-suppressed decay mode — K'^K~7t^ using a Dalitz plot tech- 
nique and find the strong phase difi'erence 5d = ^k*-k+ ~ ^K*+K- — 332° it 8° it 11° and relative 
amplitude rD = 0'K*-k+ I ^k*+k- = 0.52 ± 0.05 ± 0.04. This measurement indicates signif- 
icant destructive interference between D° — if+(if~7r°)^*- and D° — if~(if+7r'^)x*+ in the 
Dalitz plot region where these two modes overlap. This analysis uses 9.0 fb~^ of data collected at 
~ 10.58 GeV with the CLEO III detector. 
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The determination of the Cabibbo-Kobayashi-Maskawa (CKM) angle 7 (also referred to 
as 03) is important, yet challenging. Currently 7 is inferred to be 58.6°^|;9o from various 
experimental and theoretical constraints 0. Grossman, Ligeti, and Soffer |2] have proposed 
a method for a direct measurement of 7 by studying DK^, where the neutral D 

meson (D^/D^) decays to K*^K^ or K*~K^ . An important ingredient in this analysis is 
the knowledge of the relative complex amplitudes of K*^K~ and K*^K~, 

which, in the absence of CP violation, is the same as that between D° K*~K+ and 
K*^K~ . The main goal of the analysis described here is to measure the strong phase 
difference 5^ and relative amplitude between K*~K^ and K*^K~ , which 

is required for the proposed extraction of 7. We are further motivated by a recent paper of 
Rosner and Suprun [s^ that points out the sensitivity to 6d using —>■ K~^K~7i^ produced 
in e~^e~ — > iIj{3770) — > D^D^, though the analysis presented here relies on mesons 
from D*^ meson decays in e"'"e~ continuum production at ^/s ~ 10.58 GeV. This is the first 
analysis of the resonant substructures of — > K^K^tt^ and their interference. The relevant 
pubhshed individual branching ratios (BR) are BR{D^ K+K'n^) = (0.13 d= 0.04)%, 



BR{D'^ K*+K-) = (0.37±0.08)%^ BR{D^ K*' K+) = (0.20±0.11)%, and BR{D^ 
ivrO) = (0.076 ±0.005)% BSHQi. 

Three- body decays of D mesons are expected to be dominated by resonant two body 

can 



inree-Doay aecays 01 u mesons are expectea to De aommatea Dy resonant two t 
decays B H Q B and the well established Dalitz plot analysis technique Q 
be used to explore their relative amplitudes and phases. The CLEO collaboration has 
published Dalitz plot analyses for several three-body decays over the past few years 
I19L l20j and the work described here closely follows the methods developed 
in these previous analyses. 

This analysis uses an integrated luminosity of 9.0 fb~^ of e^e~ collisions at y/s ~ 
10.58 GeV provided by the Cornell Electron Storage Ring (CESR). The data were col- 
lected with the CLEO III detector B H Hi]. To suppress backgrounds and to tag the 



flavor D^{D^), the mesons are reconstructed in the decay sequence D*^ nfD^, where 
the sign of the slow pion vr+(7r7) tags the flavor of the Z}°(D°) at the time of its production. 

The detected charged particle tracks must reconstruct to within 5 cm of the interaction 
point along the beam pipe and within 5 mm perpendicular to the beam pipe (the typical 
beam spot is 300 fim in the horizontal dimension, 100 /im in the vertical dimension, and 
10 mm in the longitudinal dimension). The cosine of the angle between a track and the 
nominal beam axis must be between —0.9 and 0.9 in order to assure that the particle is 
in the flducial volume of the detector. The vr^ candidates are required to have momenta 
150 < Ptts < 500 MeV/c, and kaon candidates are required to have momenta 200 < Pk < 
5000 MeV/c. Candidate kaon tracks that have momenta greater than or equal to 500 MeV/c 
are selected based on information from the Ring Imaging Cherenkov (RICH) detector [2^ if 
at least four photons associated with the track are detected. The pattern of the Cherenkov 
photon hits in the RICH detector is flt to both a kaon and a pion hypothesis, each with its 
own likelihood Ck and Ct,. We require (— 21n£j^) — (— 21n£^) < for a kaon candidate to 
be accepted. Candidate kaon tracks without RICH information or with momentum below 
500 MeV/c are required to have speciflc energy loss in the drift chamber within 2.5 standard 
deviations of that expected for a true kaon. 

The 7r° candidates are reconstructed from all pairs of electromagnetic showers that are 
not associated with charged tracks. To reduce the number of fake 7r°s from random shower 
combinations, we require that each shower have an energy greater than 100 MeV and be in 
the barrel region of the detector. The two photon invariant mass is required to be within 
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FIG. 1: Distribution of (a) rriK+K-nO for \AM\ < 1 MeV/c^ and (b) AM for 1.84 < m^^+^-^o < 
1.89 GeV/c^ after passing all other selection criteria discussed in the text. The solid curves show 
the results of fits to the rnx+K-n" distributions, respectively. The vertical lines in (a) 

and the left-most set of vertical lines in (b) denote the signal region. The right-most set of vertical 
lines in figure (b) denote the AM sideband used for estimation of the background shape. 



2.5 standard deviations of the known tt" mass. To improve the resolution on the 7r° three- 
momentum, the 77 invariant mass is constrained to the known 7r° mass. 

We reconstruct the decay chain D*~^ ir'^D^, — > K^K~7r^ with the requirement that 
the D*^ momentum be at least as large as one-half of its maximum allowed value in order to 
suppress large combinatoric backgrounds from i?-meson decays. The candidate invariant 
mass rriK+K-nO and invariant D*~^ — mass difference AM = ^k+k--w°tt+ ~ "^k+k-xO — 
{m£)*+ —mjjo) are calculated for each candidate, where m£i*+ and m^jo are taken from Ref. j^. 
The distributions of mx+K-TvO and AM are shown in Fig. ^ We fit each of the distributions 
to the sum of two bifurcated Gaussians plus a background shape which is constant (for 
rriK+K-no) or parabohc (for AM). We find an average signal fraction of (81.8 ± 6.3 ± 2.8)%, 
where the systematic error is half of the difference between the signal fraction from the fits 
to mx+K-nO and AM. We select a signal region defined by 1.84 < mx+K-TrO < 1-89 GeV/c^ 
and \AM\ < 1 MeV/c^ which contains 735 K+K--n^ candidates. 

We expect CP violation in D decay to be negligible and assume the amplitudes for 

— s> K*^K^ and K*^K~ are equal to the amplitudes for charge-conjugated modes 

—>■ K*^K~ and K*~K~^, respectively. This allows us to double our statistics 

in a single Dalitz plot by combining flavor-tagged —>■ K^K^tt^ and —>■ K^K't:^ 
candidates and choosing the m|-_^o variable for one to be the m^+^o variable for the other 
(and vice versa). The inclusion of charge conjugate modes is implied throughout this paper. 

Figure Efa) shows the Dalitz plot distribution for the — > K^K~tt^ candidates. The 
enhanced bands perpendicular to the rnj^^^o and ^\+^o axes at an invariant mass-squared 
of ~ 0.8 GeV^/c^ correspond to i^*(892)~ and ii'*(892)"'" resonances, respectively. The 
0(1020) can be seen as a diagonal band along the upper right edge of the plot. The vector 
nature of these resonances is evident from the depleted region in the middle of each band. 
The nearly missing bottom lobe of the i^'*(892)~ band and the enhanced left lobe of the 
ir*(892)+ band show that these resonances are interfering with opposite phases with the 
S'-wave amplitude under these resonances. 



m 



We parameterize the K^K vr*^ Dalitz plot following the methodology described 

17, m, [l^ 11^. We express the amplitude for to the j*^ quasi-two-body 



state as ajC^^B) , where a,- is real and positive and B) is the Breit-Wigner amplitude for 
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FIG. 2: (a) The Dalitz plot distribution for K^K^tt^ candidates, (b)-(d) Projections onto 

of the results of Fit A showing both the fit (curve) and the 
binned data sample. The curves of Fit B projections are indistinguishable from those of Fit A. 

resonance j with spin k described in Ref. Our sign convention implies that So = 

Sk*-k+ — ^K'+K- = 0° (180°) indicates maximal destructive (constructive) interference 
between the K* amplitudes. We consider thirteen resonant components (see Table H} as 
well as a uniform non-resonant contribution. Dalitz plot analyses are only sensitive to 
relative phases and amplitudes, hence we may arbitrarily define the amplitude and phase 
for one of the two-body decay modes. The mode with the largest rate, K*^K~ , is assigned 
an amplitude aK*+K- = 1 and phase Sk*+k- = 0°. 

The efficiency for the selection requirements described above is not expected to be uniform 
across the Dalitz plot because of the momentum dependent reconstruction algorithms near 
the edge of phase space. To study these variations, we produce Monte Carlo generated 
vr+D°, D° K^K'n^ events (based on GEANT3 which uniformly populate 
the allowed phase space and pass them through our event processing algorithms. We observe 
a modest and smooth dependence of reconstruction efficiency on Dalitz plot position and fit 
this to a two dimensional cubic polynomial in (m^+^o, '^'\c+k-)- '^^^ average reconstruction 
efficiency for the decay chain D*^ tt^D^, K^K~ti^ in our signal region is found 

to be (5.8 ±0.1)%. 

Figure Q shows that the background is significant. To construct a model of the back- 
ground shape, we consider events in the data sideband 3 < AM < 10 MeV/c^ within the 
rriK+K-TTO signal region defined above, as shown in Fig. ^^b). There are 384 events in this 
selection, about three times the amount of background we estimate in the signal region. 
The background is dominated by random combinations of unrelated tracks and showers. 
Although the background includes K*^ and (p niesons combined with random tracks and/or 
showers, these events will not interfere with each other or with resonances in the signal. The 
background shape is well fitted by a two dimensional cubic polynomial in (m^+^o, ^"k+k-) 
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Resonance r 


mr (GeV/c2) 


(GeV/c2) 




0.8917 


0.0508 


0(1020) 


1.0190 


0.0043 


non-resonant 


flat 


flat 


ao(980) 


0.9910 


0.0690 


/o(1370) 


1.3500 


0.2650 


i^o(1430)± 


1.4120 


0.2940 




1.4260 


0.0985 


/o(1500) 


1.5070 


0.1090 


/2(1525) 


1.5250 


0.0730 




0.8780 


0.4990 



TABLE I: The masses and widths of resonances r considered in this analysis 



Amplitude 



Phase 



Fit Fraction 



Fit A 

SL = 18, 



6% 



K*+ 


1 (fixed) 


(fixed) 


46.1 ±3.1 


K*- 


0.52 ±0.05 ±0.04 


332 ± 8 ± 11 


12.3 ±2.2 




0.64 ±0.04 


326 ±9 


14.9 ± 1.6 


NR 


5.62 ± 0.45 


220 ±5 


36.0 ±3.7 


Fit B 








SL = 17.2% 








K*+ 


1 (fixed) 


(fixed) 


48.1 ±4.5 


K*- 


0.52 ±0.05 


313 ±9 


12.9 ±2.6 


4> 


0.65 ± 0.05 


334 ± 12 


16.1 ± 1.9 


K+ 


1.78 ± 0.43 


109 ± 17 


12.6 ±5.8 




1.60 ±0.29 


128 ± 17 


11.1 ±4.7 



TABLE II: Dalitz plot fit results. The model for Fit A includes K*^, (/), and a non-resonant 
contribution. The model for Fit B includes K*^, (p, and k^. A significance level (SL), calculated 
by the method of Ref. [2^, is shown for each fit. 

with non-interfering terms that represent K*^ and mesons. 

We use the background and efficiency parameterizations in our Dalitz plot fit to the 
data. Our results are presented in Table |n] Fit A includes the K*{892)^ and 0(1020) 
resonances plus an interfering non-resonant (NR) component and is shown in Fig. |2fb)-(d). 
For each entry in Table m the first error shown is statistical. Systematic errors are also 
shown for the K* submodes, since those are the results that ultimately contribute to the 
phase difference and relative amplitudes this analysis seeks to measure. The determination 
of these systematic errors is discussed below. 

Since it is difficult to distinguish a simple NR contribution from a broad S-wave compo- 
nent, we investigated the effect of replacing the NR component of Fit A with broad S'-wave 
— » K^n^ resonances parameterized using Breit-Wigner amplitudes 0|. The result of 
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this substitution is shown as Fit B in Table |n] Both Fit A and Fit B have good significance 
levels, and the projections of Fit A and Fit B are indistinguishable, hence we have no reason 
to prefer one fit over the other. Significance levels are calculated by the method of Ref. 

We tested other combinations of broad amplitudes as possible replacements to the simple 
non-resonant component, including one fit with Kq{1430)^ K^n^ and k,^ K^ir^ and 
another fit with a NR component combined with Kq{1430)^ . We did not find that either 
of these fits were preferable to Fit A or Fit B, although we do include these results when 
determining our model systematic error. We did not find significant evidence for any of 
the other resonances listed in Table H] A fit which included only K*^ and contributions 
(without a NR component) was significantly worse than Fit A. 

Since the choice of normalization, phase convention, and amplitude formalism may not 
always be identical for different experiments, fit fractions are reported in addition to am- 
plitudes. The fit fraction is defined as the integral of a single component (resonant or 
non-resonant) over the Dalitz plot, divided by the integral of the coherent sum of all com- 
ponents over the Dalitz plot [l5|. The sum of the fit fractions for all components will not 
necessarily be unity because of interference in the coherent sum. 

We use the full covariance matrix from Fit A and Fit B to determine the statistical errors 
on the fit fractions and to properly include the correlated components of the uncertainty on 
the amplitudes and phases. After each fit, the covariance matrix and final parameter values 
are used to generate a large number of sample parameter sets. Fit fractions are calculated as 
described above for each set of parameters, and the Gaussian widths of these distributions 
represent the statistical errors on the nominal fit fractions. 

The strong phase difference 6e> and relative amplitude r^, are defined as follows: 

where in Eq. (P) is defined as real and positive. The strong phase difference is equivalent 
to the overall phase difference due to our assumption that CP violation in D decays is 
negligible. With this definition we can simply read our nominal results from Fit A of 
Table HH 

6d = 332° ± 8° ± 11°, vd = 0.52 ± 0.05 ± 0.04. 

We consider systematic errors from experimental sources and from the decay model sep- 
arately. Contributions to the experimental systematic uncertainties arise from our models 
of the background, the efficiency, the signal fraction, and the event selection. Our general 
procedure is to change some aspect of the analysis and interpret the change in the values of 
the amplitude ratio td and phase difference 6d as an estimate of the associated systematic 
uncertainty. In Fit A, we fix the coefficients of the background parameterization to the val- 
ues found in our fit to the sideband region as described above. To estimate the systematic 
uncertainty on this background shape, we perform a fit where these coefficients are allowed 
to float constrained by the covariance matrix of the background fit. A similar method is 
used to determine the systematic uncertainty for the efficiency shape. We change selection 
criteria in the analysis to test whether our Monte Carlo simulation properly models the 
efficiency. We vary the minimum vr*^ daughter energy, the cuts on mx+K-nO and AM, the 
D*^ minimum momentum fraction, the ^^(77) — m{7r^) requirement, and the RICH and 
specific energy criteria. We allow the width of the 0(1020) to float to accommodate detector 
resolution effects. We performed partial fits of the Dalitz plot excluding regions not close 
to the i^*(892) bands, and we changed the invariant mass-squared variables in our fits from 
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{rn^+^o, '^'k+k-) ("^x-ttO' ■'^/c+tto)- The largest experimental systematic uncertainties 
are ±8° for 5d when allowing the background parameters to float, and ±0.05 for re when 
allowing the efficiency parameters to float, as described above. 

The model systematic error arises from uncertainty in the choice of resonances used to fit 
the Dalitz plot. We flt the data to many models that incorporate various combinations of 
the resonances hsted in the lower part of Table Qin addition to the K*(892)=^ and 0(1020). 
We allow our k mass and width to float in a separate flt, flnding the preferred values to be 
= (855 ± 15) MeV and = (251 ± 48) MeV. The signiflcance level for the flt where the 
K mass and width float is 16.2%. The floating k mass is consistent with Ref. jl^, but the 
floating K, width is smaller by about two standard deviations. We use our measured error 
and the error from Ref. [26j on to calculate the deviation. 

We determine the total experimental and model systematic uncertainties separately. We 
take the square root of the sample variance of the amplitudes and phases from the nominal 
result compared to the results in the series of flts described above as a measure of the 
systematic uncertainty. We flnd ±2.9°(exp.) ± 10.6°(model) for 6o and ±0.016(exp.) ± 
0.038(model) for r^j. Adding systematic errors in this way results in a model systematic 
error for Sd that is less than the difference in Sd when comparing Fit A to Fit B. We 
add the experimental and model systematic uncertainty in quadrature to obtain the total 
systematic uncertainty reported in Table |Tll 

Our systematic error is dominated by the model dependence, and the largest deviations 
from the nominal flt were observed in the series of flts where we replaced the non-resonant 
contribution with the k^. If flts including a resonance are removed from consideration, 
then the systematic errors on 6o and decrease from ±11° and ±0.04 to ±8° and ±0.03, 
respectively, and the remaining systematic uncertainty is dominated by flts including the 
iro(1430)±. 

As a cross-check, we estimate the branching ratio of K^K^tt^ from our data and 

compare it to the published value. Branching ratio measurements are not the focus of this 
analysis, so systematic errors have not been investigated. Based on our rriK+K-nO flt, we have 
a total of 627 ± 30 signal events. We may estimate the total number of D^s expected from 
continuum D*^s in our data sample, based on the integral cross-section for continuum D*^ 
production near = 10.6 GeV 0: a{e+e- D*+X) = (583 ±8 ±33 ±14) pb, where the 
fourth error stems from external branching fraction uncertainties. From this information 
and BR{D*+ vr+D") 3, we estimate BR{D'^ K+R-n'^) = (0.30 ± 0.02)%, which 
is signiflcantly higher than the previous measurement 0, Combining our flt fractions 
with known values for BR{K*^ K^tt^) and BR{(f) K^K~) we also estimate 
branching ratios of the resonant decay modes. We flnd BR{D^ — > K*^K^) = (0.38±0.04)%, 
Si?(DO ^ K*-K+) = (0.10 ± 0.02)%, and 5i?(D0 ^ ^vr^) =J0.084 ± 0.012)%. These 
branching ratios are consistent with published measurements 0, IE • 

t/-spin symmetry [3l| predicts the following for decays to a pseudoscalar meson and 
a vector meson: 

^(D° ^ TT+p-) = -A{D^ K+K*-) (2) 

and 

^(D° ^ 7r-p+) = -A{D^ -> K-K*+), (3) 
where A is the respective dimensionless invariant amplitude for each decay. Dividing Eq. (j21) 
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by Eq. ©, and assuming that A{K*+ K+n^) = A{K*- R-n^), gives: 

e*(''if*-K+-^A-+A'-). (4) 



A{D'^n+p-) aK*-K+ 



Assuming a phase convention such that a phase difference of 0° indicates maximal de- 
structive interference between p~ and p+, and assuming A{p^ 7r+7r°) = A{p^ — >■ 7r~7r°), 
we can use the recently published results of a Dalitz plot analysis of tt+tt^tt^ jl^ to 

evaluate the left hand side of Eq. (jH) : 

(0.65 ± 0.03 ±0.02)e^(356°±3°±2°) 

which may be compared to the right hand side of Eq. (@)) which comes from this analysis: 

(0.52 ± 0.05 ±0.04)e*(332°±8°±ii°)_ 

In conclusion, we have examined the resonant substructure of the decay — > K^K^tx^ 
using the Dalitz plot analysis technique. We observe resonant K*~^K~, K*~K^ , and 07r° 
contributions. We also observe a significant S-wave modeled as a n^K^ or a non-resonant 
contribution. We determine Sd = 332° ± 8° ± 11° and rn = 0.52 ± 0.05 ± 0.04. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, and the Natural Sciences and Engineering Research Council of Canada. 
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